﻿#include "LunarYear.h"
#include "ShouXingUtil.h"
#include "Solar.h"
#include "LunarMonth.h"
#include "Lunar.h"
#include <algorithm>
#include <array>
#include <cmath>
//#include <corecrt_math_defines.h>

const std::vector<std::string> LunarYear::YUAN = { "下", "上", "中" };

const std::vector<std::string> LunarYear::YUN = { "七", "八", "九", "一", "二", "三", "四", "五", "六" };

std::set<int> LunarYear::LEAP_11 = { 75, 94, 170, 265, 322, 398, 469, 553, 583, 610, 678, 735, 754, 773, 849, 887, 936, 1050, 1069, 1126, 1145, 1164, 1183, 1259, 1278, 1308, 1373, 1403, 1441, 1460, 1498, 1555, 1593, 1612, 1631, 1642, 2033, 2128, 2147, 2242, 2614, 2728, 2910, 3062, 3244, 3339, 3616, 3711, 3730, 3825, 4007, 4159, 4197, 4322, 4341, 4379, 4417, 4531, 4599, 4694, 4713, 4789, 4808, 4971, 5085, 5104, 5161, 5180, 5199, 5294, 5305, 5476, 5677, 5696, 5772, 5791, 5848, 5886, 6049, 6068, 6144, 6163, 6258, 6402, 6440, 6497, 6516, 6630, 6641, 6660, 6679, 6736, 6774, 6850, 6869, 6899, 6918, 6994, 7013, 7032, 7051, 7070, 7089, 7108, 7127, 7146, 7222, 7271, 7290, 7309, 7366, 7385, 7404, 7442, 7461, 7480, 7491, 7499, 7594, 7624, 7643, 7662, 7681, 7719, 7738, 7814, 7863, 7882, 7901, 7939, 7958, 7977, 7996, 8034, 8053, 8072, 8091, 8121, 8159, 8186, 8216, 8235, 8254, 8273, 8311, 8330, 8341, 8349, 8368, 8444, 8463, 8474, 8493, 8531, 8569, 8588, 8626, 8664, 8683, 8694, 8702, 8713, 8721, 8751, 8789, 8808, 8816, 8827, 8846, 8884, 8903, 8922, 8941, 8971, 9036, 9066, 9085, 9104, 9123, 9142, 9161, 9180, 9199, 9218, 9256, 9294, 9313, 9324, 9343, 9362, 9381, 9419, 9438, 9476, 9514, 9533, 9544, 9552, 9563, 9571, 9582, 9601, 9639, 9658, 9666, 9677, 9696, 9734, 9753, 9772, 9791, 9802, 9821, 9886, 9897, 9916, 9935, 9954, 9973, 9992 };

std::set<int> LunarYear::LEAP_12 = { 37, 56, 113, 132, 151, 189, 208, 227, 246, 284, 303, 341, 360, 379, 417, 436, 458, 477, 496, 515, 534, 572, 591, 629, 648, 667, 697, 716, 792, 811, 830, 868, 906, 925, 944, 963, 982, 1001, 1020, 1039, 1058, 1088, 1153, 1202, 1221, 1240, 1297, 1335, 1392, 1411, 1422, 1430, 1517, 1525, 1536, 1574, 3358, 3472, 3806, 3988, 4751, 4941, 5066, 5123, 5275, 5343, 5438, 5457, 5495, 5533, 5552, 5715, 5810, 5829, 5905, 5924, 6421, 6535, 6793, 6812, 6888, 6907, 7002, 7184, 7260, 7279, 7374, 7556, 7746, 7757, 7776, 7833, 7852, 7871, 7966, 8015, 8110, 8129, 8148, 8224, 8243, 8338, 8406, 8425, 8482, 8501, 8520, 8558, 8596, 8607, 8615, 8645, 8740, 8778, 8835, 8865, 8930, 8960, 8979, 8998, 9017, 9055, 9074, 9093, 9112, 9150, 9188, 9237, 9275, 9332, 9351, 9370, 9408, 9427, 9446, 9457, 9465, 9495, 9560, 9590, 9628, 9647, 9685, 9715, 9742, 9780, 9810, 9818, 9829, 9848, 9867, 9905, 9924, 9943, 9962, 10000 };

std::vector<int> LunarYear::YMC = { 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 };

LunarYear::LunarYear(int lunarYear)
{
	Year = lunarYear;
	int offset = lunarYear - 4;
	int yearGanIndex = offset % 10;
	int yearZhiIndex = offset % 12;
	if (yearGanIndex < 0)
	{
		yearGanIndex += 10;
	}
	if (yearZhiIndex < 0)
	{
		yearZhiIndex += 12;
	}
	GanIndex = yearGanIndex;
	ZhiIndex = yearZhiIndex;
	Compute();
}

LunarYear* LunarYear::cacheYear = nullptr;


const LunarMonth* LunarYear::GetMonth(int lunarMonth) const
{
	auto it = std::find_if(Months.begin(), Months.end(),
		[this, lunarMonth](const LunarMonth& m) {
			return m.Year == this->Year && m.Month == lunarMonth;
		});
	return (it != Months.end()) ? &(*it) : nullptr;
}


void LunarYear::Compute()
{

	// 1. 初始化数据结构
	std::array<double, 27> jq = { 0 }; // 27个节气
	std::array<double, 16> hs = { 0 }; // 16个合朔时刻
	std::array<int, 15> dayCounts = { 0 };
	std::array<int, 15> months = { 0 };

	const int currentYear = Year;
	// 2. 基础天文计算
	const double jd = std::floor((currentYear - 2000) * 365.2422 + 180);
	double w = std::floor((jd - 355 + 183) / 365.2422) * 365.2422 + 355;
	if (ShouXingUtil::CalcQi(w) > jd) {
		w -= 365.2422;
	}
	// 25个节气时刻(北京时间)，从冬至开始到下一个冬至以后
	for (auto i = 0; i < 26; i++)
	{
		jq[i] = ShouXingUtil::CalcQi(w + 15.2184 * i);
	}

	// 4. 节气精确计算（从大雪到立春）
	JieQiJulianDays.clear();
	for (int i = 0; i < Lunar::JIE_QI_IN_USE.size(); ++i) {
		double jdVal;
		if (i == 0) {
			jdVal = ShouXingUtil::QiAccurate2(jq[0] - 15.2184);
		}
		else if (i <= 26) {
			jdVal = ShouXingUtil::QiAccurate2(jq[i - 1]);
		}
		else {
			jdVal = ShouXingUtil::QiAccurate2(jq[25] + 15.2184 * (i - 26));
		}
		JieQiJulianDays.push_back(jdVal + Solar::J2000);
	}
	// 5. 合朔计算
	w = ShouXingUtil::CalcShuo(jq[0]);
	if (w > jq[0]) {
		w -= 29.53;
	}

	for (int i = 0; i < 16; ++i) {
		hs[i] = ShouXingUtil::CalcShuo(w + 29.5306 * i);
	}
	// 6. 月天数计算
	for (int i = 0; i < 15; ++i) {
		dayCounts[i] = static_cast<int>(hs[i + 1] - hs[i]);
		months[i] = i;
	}

	// 7. 闰月判定
	int prevYear = currentYear - 1;
	int leapIndex = 16; // 默认无闰月

	if (LEAP_11.find(currentYear) != LEAP_11.end()) {
		leapIndex = 13;
	}
	else if (LEAP_12.find(currentYear) != LEAP_12.end()) {
		leapIndex = 14;
	}
	else if (hs[13] <= jq[24]) {
		int i = 1;
		while (i < 13 && hs[i + 1] > jq[2 * i]) {
			++i;
		}
		leapIndex = i;
	}

	// 8. 月份调整
	for (int i = leapIndex; i < 15; ++i) {
		months[i] = months[i] - 1;
	}

	// 9. 农历月构建
	Months.clear();
	int fm = -1;
	int index = -1;
	int y = prevYear;

	for (int i = 0; i < 15; ++i) {
		const double dm = hs[i] + Solar::J2000;
		int v2 = months[i];
		int mc = YMC[v2 % 12];

		// 特殊年份修正（需补充实际天文历法规则）
		if (dm >= 1724360 && dm < 1729794) {
			mc = YMC[(v2 + 1) % 12];
		}
		else if (dm >= 1807724 && dm < 1808699) {
			mc = YMC[(v2 + 1) % 12];
		}
		else if (std::abs(dm - 1729794) < 1e-9 ||
			std::abs(dm - 1808699) < 1e-9) {
			mc = 12;
		}

		if (fm == -1) {
			fm = mc;
			index = mc;
		}

		if (mc < fm) {
			y += 1;
			index = 1;
		}
		fm = mc;

		// 闰月标记
		if (i == leapIndex) {
			mc = -mc;
		}
		else if (std::abs(dm - 1729794) < 1e-9 ||
			std::abs(dm - 1808699) < 1e-9) {
			mc = -11;
		}

		Months.emplace_back(y, mc, dayCounts[i], dm, index);
		++index;
	}
}






const LunarYear& LunarYear::FromYear(int lunarYear)
{
	LunarYear* y;
	if (cacheYear == nullptr || cacheYear->Year != lunarYear)
	{
		delete cacheYear;
		cacheYear = nullptr;
		cacheYear = new  LunarYear(lunarYear);
		y = cacheYear;
	}
	else
	{
		y = cacheYear;
	}
	return *y;

	//return new  LunarYear(lunarYear);
}


















